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1. Introduction 



The infrared behaviour of the QCD Green's functions contains essential information about the 
realisation of confinement in the covariant formulation of QCD, especially in Landau gauge [jjj]. It 
has been focus of intensive research over the past ten or so years and was studied with a variety 
of methods. These include different continuum-based functional methods ^J, ||, ||] as well as 
Monte Carlo (MC) simulations. In particular, the former have predicted that, in Landau gauge, 
the deep-infrared behaviour of the momentum-space gluon and ghost dressing functions, Z and G, 
should be given by [^] 

Z{p 2 ) ~ (//Aqcd) 2 ^ and G(p 2 ) ~ {p 2 /^qcd) Kg forp 2 ^0, (1.1) 

which are both determined by a unique critical infrared exponent Kz = Kq = K £ (0.5, 1). Under a 

mild regularity assumption on the ghost-gluon vertex [Qj, the value of this exponent is furthermore 

obtained as M, |] , 

k = (93- Vl20T)/98 ss 0.595. (1.2) 

The prediction of this so-called conformal infrared behaviour has been challenged by a number 
lattice simulations. In fact, with the notable exception of pure SU(2) lattice gauge theory in 2 
dimensions |jg], the proposed infrared behaviour ( |l.l| ), or even its onset, has not conclusively been 
confirmed in lattice simulations so far. In contrast, numerical evidence is provided that the dressing 
functions are consistent with a transverse gluon mass, Z ~ p 2 /M 2 and G = const., for p 2 — > 0. 

Finite-volume effects were soon identified as the prime suspects responsible for the disagree- 
ment between continuum and lattice results [0]. In fact, in order to be able to see an onset of an at 
least approximate conformal infrared behaviour in a finite volume of extend L, a reasonably large 
range of momenta satisfying 

n/L < p < Aqcd (1.3) 

needs to be accessible, and it was concluded in [|7|] that lattice sizes of about 15 fm are needed before 
even an onset of the leading infrared behaviour can be observed. Not only this has triggered a run 
for the biggest-ever lattice sizes used in simulations (see, e.g., [||). All those simulations, however, 
show very little if no tendency to follow the finite-size corrections to the predicted asymptotic 
behaviour, even though their volumes reached up to about 20 fm in size. 

In this study, we focus on the gluon and ghost propagator in the strong coupling limit, j8 — ► 0, 
of pure SU (2) lattice Landau gauge. This unphysical limit, which can be interpreted as the formal 
limit Aqcd — ► °°> allows us to assess whether the predicted conformal behaviour can be seen for 



the larger lattice momenta p, after the upper bound in (|L3|) has been removed, in a range where 
the dynamics due to the gauge action would otherwise dominate and cover it up completely. Our 
complete analysis can be found in [0]. 



2. Gluon and ghost dressing functions in the strong coupling limit 

We simulate SU(2) gauge theory in the strong coupling limit by generating random link con- 
figurations {£/}. These are sets of SU(2) gauge links, U xfl = u®l + io a u"^, equally distributed 
over (u P,u) X fi G S 3 . Those configurations are fixed to the (standard lattice) Landau gauge (SLG) 
using an over-relaxation algorithm that iteratively minimises the SU(2) gauge functional of SLG, 
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Figure 1: The gluon propagator (left) and its dressing function (right) vs. lattice momentum a 2 q 2 for dif- 
ferent lattice sizes in the strong coupling limit. Straight lines (right) represent fits of the 56 4 data to its 
asymptotic behaviour at low (solid) and high (dashed) momenta (see |]]). 



V £ 7[g]=4£(l-iTrf/| i ) . (2.1) 

To satisfy Landau gauge with sufficient accuracy the over-relaxation algorithm is iterated until 
the stopping criterion max^Tr [(V*Af^)(V*Af^)] < 10~ 13 is met at every lattice site. Here Vj^ 
denotes the lattice backward derivative and Af „ is the lattice gluon field of SLG, given by 

in terms of the gauge-transformed link UL. Gluon and ghost propagators are then calculated in 
momentum space employing standard techniques. For further details refer to [QJ. 

The data for the gluon propagator and its dressing function is shown in Fig. [I] versus the 
lattice momentum a 2 q 2 . The propagator is seen to increase with momentum, while it plateaus 
at low momenta. Perhaps unexpectedly, however, this happens irrespective of the lattice size at 
around a 2 q 2 w 1. It is therefore not a finite- volume effect, and the observed mass behaves as 
D 1 (0) oc M 2 oc l /a 2 in the strong-coupling limit with no significant dependence on L. 

The asymptotic behaviour at either end is also seen in the corresponding dressing function Z. 
We compare this function to fits obtained from the 56 4 data in the right panel of Fig. [[]: Writing 
x = a 2 q 2 , we use f(x) = cx 2Kz as the fitting function for x £ [2, 14] and f(x) = dx for small momenta 
x < 0.1 to fit c, d and the exponent Kz- The data at large momenta is well described by the power 
law ( JTTTt ) with an exponent Kz = 0.563(1), while it approaches the massive behaviour at small x. 

To assess the model dependence of these results we also used different fit models and a broad 
range of fitting windows (see [^] for details). Monitoring the corresponding x 2 / ndf values we find 

that , , „ 

D I (x)=cx 2Kz - [ and D n (x) = c{d + x) 2Kz - 1 (2.3) 

provide stable fits with similar values for Kz and very little dependence on the lattice size. While 
Dj(x) performs slightly better at large momenta, D n {x) provides good stable fits of the gluon 
propagator over the whole momentum range. 
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The Kz values extracted from dif- 
ferent fit models and lattice sizes with 
fixed window x G [2,14] are shown in 
Fig. ||. The observed dependencies on 
either one are rather small. There is a 
general trend for Kz to slightly increase 
with a/L (the dark grey band in Fig. ^ 
though this is within the systematic un- 
certainty due to the fit model. 

Similar fits were performed to ex- 
tract the exponent k g from the ghost 
dressing function G. These fits are less 
robust with a more pronounced model 
dependence (light grey band in Fig. ^). 
This is mainly due to the wider transi- 
tion region, from G = const, at small x 

to G ~ x~ K ° at large x, which is under less control here. The exponent can nevertheless be estimated 
as Kq = 0.60(7). The results are consistent with the scaling relation K z = K G . 
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Figure 2: K versus a/L for the gluon and ghost propagators. 



3. Comparing different lattice definitions of gauge fields 

It is interesting to compare the above presented SLG data to that of the modified lattice Landau 
gauge (MLG) introduced in [ |I0| ] which is based on stereographic projection to define lattice gauge 
fields. When comparing MLG to the ever popular SLG, there is no advantage that the SLG has over 
the MLG. A promising particular feature of the MLG on the other hand is that it provides a way 
to perform gauge-fixed MC simulations sampling all Gribov copies of either sign (of the Faddeev- 
Popov determinant) in the spirit of BRST. This feature will be explored in a forthcoming study. 
Here we simply use the MLG for comparison in the standard way, i.e., we gauge-fix configurations 
via minimisation of the MLG functional, which in SU(2) is of the form [ |To| ] 

%fe] = -8£lnQ + I T r^ A ) . (3.1) 



The MLG functional differs from that of SLG (2. 1 ) and so do the lattice gluon fields of MLG, 

1 ~+\ . ~ 2U, 

1 1 

2 11 ^xn 



A., = -^-^) where u^*^. (3.2) 



The lattice Landau gauge condition F x = V* A XjU = is formally unchanged in the stereographically 
projected gauge fields A, in particular, we have F X (A) = F X (A). The corresponding Faddeev-Popov 



operator is given explicitly in [10]. Both lattice definitions of Landau gauge have the same con- 
tinuum limit, and any differences between MLG and SLG data at finite lattice spacings are lattice 
artifacts. It is also worth mentioning that gauge configurations fixed to MLG do not satisfy the 
gauge condition of SLG and vice versa. Nonetheless, exact transversality, i.e., q^tyA^k) = 0, is 
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satisfied at finite lattice spacing a for both of them equally, if the momenta (k) are defined as 
aq fl (k) :=2sm(nk^/N^) with integer valued k^ G (—N^/2,N^/2] in the usual way. 

The data for the gluon propagator of SLG (red filled diamonds) is compared to that of MLG 
(blue filled circles) in Fig. ||. There we also show data for the gluon propagator where either 

aA$ = u%u a m a a , or aA% = <j>^G a /2 from U xfl = exp{^ a a /2} (3.3) 

are used to define lattice gluon fields based on the adjoint representation, A ad} (black open dia- 



monds), and thus blind to the centre [11], or on the tangent space at the identity A ln (green crosses). 
In these two cases, A 3dl and A ln , for the purpose of a qualitative comparison, we simply use the 
gauge configurations of the SLG to calculate the gluon propagator. Especially for A ln this implies, 
however, that the condition q^(k)A^(k) = is satisfied at best approximately and nowhere near the 
precision of the other two (SLG and MLG). The residual uncertainty due to other possible tensor 
structures then causes the somewhat larger errors for this definition as seen in Fig. ||. 



First, we fit the data from all four definitions to D u (x) in ( |2.3| ) which provides the best overall 
description in the full momentum range. In order to demonstrate how the other definitions compare 
to the SLG, we keep its value for the exponent fixed when fitting the other data, i.e., Kz = 0.562 as 
obtained for L/a = 32 in SLG is used in all fits. Relative to the scaling branch <=c x lKz ~ x for large 
x = a 2 q 2 we then observe a strong definition dependence in the (transverse) gluon mass term at 
small x (see the left panel in Fig. [|). The relative weight of the two asymptotic branches, scaling 
at large x and massive at small, is clearly discretisation dependent and can not be compensated by 
finite renormalisations. A first indication that the massive branch might indeed be the ambiguous 
one is the observed M« l/a. This is consistent with the fact that the definitions of gauge fields 
on the lattice, which agree at leading order, all differ at order a 2 , and so do their corresponding 
Jacobian factors which leads to lattice mass counter-terms of different strengths. The observation 




Figure 3: The gluon propagator in the strong-coupling limit over x — a 2 q 2 for the various definitions of 
gauge fields mentioned in the text. All data from 32 4 lattices and normalised to the scaling branch. All fits 
were performed with fixed Kz — 0.562 as obtained for the SLG (with L/a = 32). On the left, the results 
are shown for Du (x) in ( 2.3 1. On the right we use Dm =d + cx 2Kz ~ 1 to demonsttate how the data from all 
definitions collapse to a unique curve °< x 2Kz ~ l (solid, red) when the mass term is subtracted. 
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that these differences matter here explicitly demonstrates the breakdown of the lattice Slavnov- 
Taylor identities in minimal lattice Landau gauge in the non-perturbative domain. 

To assess whether this ambiguity has an influence on the exponent Kz, we have also used fits 
of the form D in = d + cx 2Kz ~\ again with the fixed K z = 0.562 from the 32 4 SLG data. This form 
leads to somewhat larger # 2 /ndf that arise from the transition region around x = a 2 q 2 ~ 1 which 
is described a bit better by Djj. If we then subtract the constant d, however, the normalised data 
of all four definitions nicely collapse onto a unique curve oc x 2Kz ~ x as seen in the right panel of 
Fig. 0. The data is fully consistent with a unique exponent Kz of around the SLG value which has 
an infinite volume extrapolation of 0.57(3). 

The strong-coupling ghost propagators of SLG and MLG are presented in [[J. The best global 
fits are then obtained for the form G~ x (x) = d + cx^. They are again consistent with a unique 
scaling exponent Kq Kz and deviate from one another in the relative strength of the constant term 
at small x = a q 2 . Instead of discussing the strong-coupling ghost propagator in more detail here 
we turn directly to the product of dressing functions defining a running coupling. 



4. Coupling constant 



The predicted infrared scaling ( |1.1| ) with Kz 
pling defined by [0] 2 



Kq immediately implies that the running cou- 

(4.1) 



fjZ(/)GV 
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a s {p) 

approaches an infrared fixed point, a s — > 
a c for p 2 — > 0. Standard continuum con- 
ventions of course need rescaling g 2 Z — > Z 
when comparing to lattice definitions such 
as (22). The predicted conformal scal- 
ing in the strong-coupling limit, with Z = 
c z {a 2 q 2 ) 2K and G 1 = c G (a 2 q 2 ) K , would 
therefore imply that the coupling ( |4.1[ ) is 
indeed constant with a s = a c = cz/{4ftc 2 G ). 
Note that its value is thus determined pre- 
cisely by those multiplicative constants in 
the propagators that have been ignored in 
the analysis up to here. They have to be 
extracted from the bare lattice data with- 
out rescaling or renormalisation. In com- 
plete agreement with the general observa- 
tion of conformal scaling in the strong- 
coupling limit at large momenta the product of the gluon and ghost dressing functions levels 
at a c » 4 for large a 2 q 2 . As expected for an exponent K slightly smaller than the value in (1.2) this 
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Figure 4: The strong-coupling limit of 0C S from (1.1) in 
minimal lattice Landau gauge for SLG and MLG. The dot- 
ted line marks a c max w 4.46, the maximal value of the criti- 
cal coupling for SU (2) from the continuum prediction [ph. 



is just below the upper bound a f max ss 4.46 for SU (2), see [Bt]. It is quite compelling that also this 
result is nearly independent of the gauge-field definition, i.e., almost identical for SLG and MLG, 
see Fig. |] Predominantly driven by the ghost propagator, the violations to this conformal infrared 
scaling set in as soon as the ambiguity in the definition of minimal lattice Landau gauge does. 
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5. Conclusions 

We have studied gluon and ghost propagators of pure SU(2) minimal lattice Landau gauge in 
the strong coupling limit. This unphysical limit probes the gauge field measure of the minimal 
lattice Landau gauge for there is no contribution from the Yang-Mills (plaquette) action. 

As expected for a formal limit Aqcd — > 00 , it is then observed that the propagators show a 
conformal scaling behaviour for large lattice momenta, a 2 q 2 S> 1 . Finite-size effects are small 
and the combined gluon and ghost data is consistent with an L/a — > °° extrapolation of a critical 
exponent K = 0.57(3). This scaling branch at large a 2 q 2 furthermore leads to a critical coupling 
of a c » 4 which is just below the predicted maximum a™ ax 4.46 for SU (2). These results show 
very little if no significant dependence on the lattice definition of gauge fields and measure. 

Another unambiguous result is the emergence of a transverse gluon mass M« l/a in the 
strong-coupling limit of minimal lattice Landau gauge. Both propagators show this massive be- 
haviour at small momenta corresponding to Z ~ q 2 /M 2 and G = const, for a 2 q 2 <C 1. This massive 
low momentum branch of the data, however, depends strongly on which lattice definition is being 
used for the gauge fields and their measure. This is typical for a mass counter-term on the lattice 
and demonstrates the breakdown of lattice Slavnov-Taylor identities (STIs) and BRST symmetry 
in minimal lattice Landau gauge beyond perturbation theory. It is still possible that this ambiguity 
disappears in the continuum limit, eventually. But because it is a combination of ultraviolet (mass 
counter-term) and infrared (breakdown of STIs) effects, this might take very fine lattice spacings in 
combination with very large volumes and therefore who-knows-how big lattices to verify explicitly. 

This research was supported by the Australian Research Council and by eResearch South Australia. 
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